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We study skew-orthogonal polynomials with respect to the weight function exp[— 2V^(a;)], with 
V{x) = ^if^iiuK / K)x^ , U2d > 0, d > 0. A finite subsequence of such skew-orthogonal polynomials 
arising in the study of Orthogonal and Symplectic ensembles of random matrices, satisfy a system 
of differential-difference-deformation equation. The vectors formed by such subsequence heis the 
rank equal to the degree of the potential in the quaternion sense. These solutions satisfy certain 
compatibility condition and hence admit a simultaneous fundamental system of solutions. 
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, 1. INTRODUCTION 

■ The concept of 'Universality' in random matrix theory and its various appUcations in real physical s yste ms have 
' attracted both mathematicians and physicists in the last few decades P, 0, H, 0, 0, 0, H, M, flllll^. fisl [13]. 
^ From mathematical point of view, the study of 'Universality' in the energy level correlations of random matrices 
[ require a good understanding of the asymptotic behavior of certain families of polynomials. For example Unitary 
I ensembles of random matrices require an understanding of the large n behavior of orthogonal (the one-matrix model) 
' and bi-orthogonal (the two matrix model) polynomials, while Orthogonal and Symplectic ensembles require that of 
, the skew-orthogonal polynomials. The rich literature available on orthogonal polynomials [1, 0, H, H, [13 Hlli [13, [Hj 
\^ • [13, m, [H, uM ^^'^ bi-orthogonal polynomials [l^, [13, [HI [H, [11] have contributed a lot in our understanding of the 
I Unitary ensembles. Our aim is to develop the theory of skew-orthogonal polynomials [1, [3, [E [13 [H) [H, [13 so that 
we can have further insight into the one-matrix model for Orthogonal and Symplectic ensembles of random matrices. 

In this direction, previous experience with orthogonal and bi-orthogonal polynomials arising in the two-matrix model 
makes us believe that perhaps the most logical and rigorous way to study the asymptotic properties of these skew- 
' orthogonal polynomials is to do the following: (1). First, we must look for a finite subsequence of skew-orthonormal 
, vectors which will satisfy a mutually compatible system of differential-difference-deformation equation. This requires 
' one to fix the 'size' or rank of this finite sub-sequence. The Generalized Christoffel-Darboux formula helps in this 
' regard as it gives an indication of the number of terms near the so-called 'Fermi-level' which are actually required to 
study this system. (2). The next step is to look for a mutually compatible system of differential-difference-deformation 
equation satisfied by the finite subsequence of these vectors and hence find the so-called 'fundamental solution'. (3). 
Finally, to formulate a quaternion matrix Riemann-Hilbert problem and by applying the steepest descent method, 
' obtain the Plancherel-Rotach-type [IB] formula for skew-orthogonal polynomials. 

Reference fsl] achieved the first goal. In this paper, we try to address the second issue of the existence of the 
'fundamental solution' of a system of differential-difference-deformation equation. One of the by-products is that 
we observe certain duality property between the skew-orthonormal vectors of Orthogonal and Symplectic ensembles. 
(In fact, this justifies further the use of skew-orthogonal polynomials to ordinary orthogonal polynomials in studying 
these ensembles [i,[l[28, 29, 30, 31, 32, 33, 34, 35, 36].) The third part, which is also the most crucial one, is not yet 
fully understood. 

We shall consider an ensemble of 2N dimensional matrices H with probability distribution 

Pp,N{H)dH = -i_expl-[2TrF(iI)]]di?, (1.1) 
where the parameter /3 = 1 and 4 corresponds to the joint probability density for Orthogonal and Symplectic ensembles 
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of random matrices. The "potential" V{x) is a polynomial of degree 2d with positive leading coefRcient: 

2d 

"2d>0, d>0, (1.2) 

A'=l 

where is called the deformation parameter. Zp^ is the so called 'partition function': 



ZpN ■■= / exp[-[2Trl-(i/)]]di/ - iV! [] gf\ (1.3) 



where M2^^ is a set of all 2N x 2N real symmetric (/? = 1) and quaternion real self dual (/3 — 4) matrices. dH 
is the standard Haar measure, gj^'' is the skew-normalization constant for polynomials related to Orthogonal and 
Symplectic ensembles (2^ . 

Remark on notation. Before entering the details of calculation, we must mention that throughout this paper, we 
have followed to a great extent notations used in [13, [H, HO) HH ■ Apart from having the advantage of using a well 
established and compact notation, this strategy will also highlight the striking similarity (as well as the difference) 
between skew-orthogonal polynomials and their bi-orthogonal counterpart. 

B. Skew-Orthogonal Polynomials: relevance in Orthogonal and Symplectic ensembles 
Definition: 

For Orthogonal and Symplectic ensembles of random matrices, we define semi-infinite vectors: 

$('5)(x) = ($[,'')*(x)...$(f)*(a;)...)*, $('5)(a;) = -$('5)*(x)Z, (1.4) 

^(/3)(a;) = (*[,'^)*(x)...*(f)*(a;)...)*, $(«(a;) = -*('^)*(x)Z, (1.5) 

where 

*W(x) = $f )(x), ¥„'\x) = ^cf>(i)(y)e(:. - y)dy. (1.6) 
Each entry in these semi- infinite vectors is a (2 x 1) matrix: 



where 



is the weighted skew-normalized polynomial (often called the quasi-polynomial) and 



is monic skew-orthogonal polynomial of order n 



-1 

is semi-infinite anti-symmetric block-diagonal matrix with = —1 and 



(1.7) 



0(f) (:e) = -^nlt\x)exp[-V{x)], (1.8) 



4^)(x) = ^ci-'^V, c("'« = l, neN, (1.9) 

fe=0 



Z ^ [ ^^l ]+...+ (1.10) 



e{r) = n (1.11) 
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is the step function. These vectors satisfy skew-orthonormaUty relation: 



1 

nrn | g ^ 



n, TO e N. 



(1.12) 



Remark This definition difi^ers from that of d i, for /3 = 4 by a factor of 2, which is incorporated in the 
normahzation constant. 

These skew-orthonormal vectors form the necessary constituents of the kernel functions required to study different 
statistical properties of Orthogonal ensembles and Symplectic ensembles of random matrices. For example, the two 
point correlation function can be written in terms of the 2x2 kernel function [26. . .27( | : 



^2Ni^j y) 



where d is the kronecker delta and 



(1.13) 



Si^j^ {x,y)^m (x) n (y) , Di%,y) = (x) [] ^^^'^ iv) , 



N 



N 



(1.14) 



i^/J {x, y) = {x) n (y), (2/, x) = m (x) n 'f^^'' (y) 



N 



N 



The level density is given by 



The matrix 



p^^Hx,x)^S^/Jix,x) 



n = diag(l,...,l,0,...,0) 



N 



N 



(1.15) 



(1.16) 



(1.17) 



is formed by (2 x 2) unit matrices (i.e. a unit matrix of size 2N in real space). In general, this will be used to 
truncate semi-infinite matrices and in the above case, the semi-infinite vectors of the finite sum in Eas. (|1.14p and 

(nni). 



The Generalised Christoffel Darboux sum 



Here, we present a summary of the main results of Ref. where we studied the kernel function S2^(x,y). We 

expand (x^^'^-' (x)), {^^^\x))' and {x^'^^\x))' in terms of $'^''^(2;) (and hence introduce the semi-infinite matrices 
Q(/3)^ p(/3) and R'^l^^ respectively): 

^ Q(^)<^>W{x), (1.18) 

*(4)(x) = P(4)$(4)(a;), x'i>^^\x)^R^^^<i>^^\x), (1.19) 

2$(i)(x) = P(i)*(i)(a;), 2x<i>'-^\x) = R''^'^¥^\x), (1.20) 

where Eq. (|1.20p is obtained by multiplying the above expansion by e(y — x) and integrating by parts. Here 

=p(4)g(4)^ ^(1) ^Q(l)p(l). (1.21) 

The matrices P'-'^\ Q'-'^^ and R'-^'^ are quaternion matrices. For a nice introduction to the subject, the reader is 
referred to the book by Prof. M. L. Mehta [13 . 
Defining the quaternion-matrix 



i?''^(x) :=i?W_3;p(4)^ 



,p(i) 



(1.22) 
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the Generalized ChristofFel Darboux formula for the Symplectic ensembles {(3 = 4) can be written as: 



N 



(1.23) 



For Orthogonal ensembles {(3 = 1) the Generalized Christoffel Darboux is given by: 



N J 



¥'Hy). 



(1.24) 



The Generalized Christoffel-Darboux matrix 
the form 



in terms of the elements of the quaternion matrix takes 



N 



-R 



\ 



■2N,2N-2d 






{x) 



-R 



'2N,2N-1 



ix) 



^2N-2d.2N\^> U 



^2N-1,2N\^) ■■■ n 

















27V-l,2Af+2(i-l' 









x) 



,iV > 41.25) 



Here, we must point out that there is a small difference with the notation used in the second part of Ref. ^ , where 
we used a matrix of size 2N + 2 to prove the 'Universality' in the Gaussian case. Also the definition of i?'-'^' and P'^^ 
differs by a factor of 2. 



D. Difference-differential-deformation equations 



A quick glance at Eq. ()1.25|) reveals that the relevant vectors contributing to the correlation function are 
'?^(2A'+2d-i)(^)' This prompts us to define a finite subsequence (or window: W ^ {N — d, . . . ,N + d~ 
1},N > d) of skew-orthogonal vectors of size 2d (in the quaternionic sense), ^^^\x) and "^^\x): 



such that 



(/3) 



\ <P2Af+2d-l 
( 4'L2di^) \ 



V '^2N+2d-li^) J 



(/3) 
N+d- 



^l(^) / 



V <id-ii^) I 



N > d. 



N> d. 



(1.26) 



(1.27) 



v.? e N, 
\/j e N. 



(1.28) 
(1.29) 



The rank of the window is equal to the degree of the potential V{x) in the quaternionic space and twice that of the 
degree of V{x) in the real space. 

The recursion relations connecting the finite subsequence (or window) with the upper or lower one is through the 
ladder operator 



$[^)+i(a;) = A^^^(x)$l^^(x), $l^^_,(x) = iA'^Zii^)r''i>^^>ix) 



l(4)/ 



,(4) 



(4) 



l(4) 



-lfl>(4). 



(1.30) 



,(1) 



(1) 



,(1)/ 



and 



s(i) 



(1) 



s(i)/ 
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(1.31) 

(1.32) 
(1.33) 



The vectors ^'^^ {x) and (x) also satisfy a system of ODEs: 



— ^^^hx) 
dx^'^^''' 


= d!-^\x)^'^\x), 




= D«(^)^1^^W, 






^ D'^^\x)<i>'^\x), 




^D^^\x)¥^\x). 


(1.34) 



Under an infinitesimal change of the deformation parameter uk (the coefficients of the polynomial potential 
these vectors satisfy a system of PDE's given by 

^^'^{x) = ut'\x)^^^{x), I ^<^(.) = f/^^^^(x)<)(.), 



V{x)), 



' ^'^{x) = U-''\x)^^^{x), 



duK 



du 



K 



^\:^{x)=u^^'\x)^'^{x). 



(1.35) 



Thus the finite subsequence of skew-orthogonal vectors satisfy a system of difference-differential-deformation equation. 

Remark. The two pairs of matrices ^j^^(x) and {x) (similarly {d'^^{x) and D^^{x)), U!k'^\x) and Uj^^^\x) 
(similarly {U^^^\x) and U^^^\x)) and A^^\x) and A^^\x) (similarly B^^\x) and B^\x)) are dual in the sense 

that they remain invariant under an interchange of ^'^^^(a;) ^'^^\x) and Q^*^ Q^^^ (and similarly for '^^^\x) and 
$(4)(^)). 



E. Compatibility 

The existence of recursion relations, differential equations and deformation equations for vectors arising in the 
Orthogonal ensembles and Symplectic ensembles of random matrices can be viewed as just a projection of the semi- 
infinite functions ^^^\x) and ^^^\x) onto the finite "window" ^^{x) and ^^^^^(a;) respectively. However, we may also 
consider these equations as defining an ovcrdctermined system of finite difference-differential-dcformation equations 
of the vector functions and see that these systems are compatible. This leads to the existence of a fundamental matrix 
solution, denoted by $j^j(a;) and 4'j^j(a;) where all the column vectors satisfy the above difference-differential- 
deformation equations simultaneously. 

The compatibility of the deformation and difference equation with the differential equation imply that the general- 
ized monodromy of the operator — D^j^\x)j and — D''^\xfj is invariant under uk deformations and shifts 
in N. 



Outline of the article 

-Section 2 deals with the properties of different finite band matrices related to skew-orthogonal polynomials. 
-In section 3, we study the system of PDE's arising from the infinitesimal change of the deformation parameter Uk- 
-We derive the difference relations satisfied by the finite subsequence of skew-orthogonal vectors in section 4. 
-In section 5, we derive the folding function which is used to project any given vector (or ^l^^) onto its finite 
subsequence or window. 

-In section 6, we obtain the folded deformation matrix, using the results of section 5. 

- The differential equation for skew-orthonormal vectors, using the results obtained in Section 6, is derived in section 

7. 

-In section 8, we discuss the existence of the Cauchy-like transforms of the skew-orthogonal vectors of order n as 

the other solutions to the diffcrcntial-deformation-difference equations, for fairly large n. 

-We prove compatibility conditions for these difference-deformation-differential system in section 9. 



6 



2. RECURSION RELATIONS AND FINITE BAND MATRICES. 

In this section, we will study in detail the different recursion relations satisfied by the skew-orthonormal vectors. 
Here, we must note that unlike the orthogonal polynomials, skew-orthonormal vectors do not satisfy a three-term 
recursion relation and hence do not give birth to the tri-diagonal Jacobi matrix. 

Differentiating Eas. pTT9)) - pr20| . we get 



d_ 

dx 



$(/3)(a;)=p(«$(«(x), 



dx 



^-(/3)(a;) =p('3)*(«(a;), /3=I,4, 



^_i$(4)(2.) ^ p(4)$(4)(^)^ ^_1^(4)(2.) ^ r^(4) _ ^(4)(^)^ 

dx dx I i 



dx dx 
where ^'^^\x) and '^^^\x) are semi-infinite column vectors defined in Eas. (|1.4ll.5p . This is equivalent to saying 



(2.1) 
(2.2) 
(2.3) 



[i?(4) _ ^p(4)]$(4)(^) ^ 0, - a;P(^) - l]*(4)(a;) - 0, 

These finite band matrices satisfy the following commutation relations: 



(2.4) 
(2.5) 



Here, each entry is a 2 x 2 quaternion of the form: 



[P(«,P(^)]=P(«. 



^2n+l,2m ^2n+l,2m+l 



Using ('(/'n'\ V'm ^ ) and (a;i/'n'\ V"™^ ) for = 4, and replacing %l)^'^\x) by (/)(^)(a;) for /3 = 1, we get 



p(/5) ^ _p(/3)^^ 



p(/3) ^ 



where dual of a matrix A is defined as [27 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



Remark. The condition of anti-self duality imposes a much lighter restriction on a diagonal quaternion than 
antisymmetry condition on a diagonal matrix element. For example, it leaves the off-diagonal entry of the diagonal 
quaternion arbitrary. This is the reason why the odd skew-orthogonal polynomials are arbitrary upto the addition of 
a lower even polynomial. 

However starting with {x4>\i'' , i/'m ^ ) and using (|2.6[) , we get 



(2.10) 



As pointed out in the beginning of the section, this essentially means that like the orthogonal polynomials, we do 
not have a tri-diagonal Jacobi-matrix for skew-orthogonal polynomials. It is this relation that causes the Generalized 
Christoffel-Darboux sum to have a local behavior (i.e. the rank has dependence on the measure or weight function). 
In a sense, this is a major setback in our hope of defining a 2 x 2 matrix Riemann-Hilbert problem for skew-orthogonal 
vectors, similar to that of the orthogonal polynomials. Also from the dependence of P'^^'^ and R'^^^ on 2d (the degree 
of the polynomial V{x))^ we can conclude that the 'size' of our q-matrix Riemann-Hilbert problem will depend on d. 
From 1^, we also have 



F^'^) + [V'iQ'^^^ - lower, R^l^'> + [Q'^^W{Q'^^^ 



: lowers 



(2.11) 



where 'lower' and 'lower+' denotes a strictly lower triangular matrix and a lower triangular matrix with the principal 
diagonal respectively. 
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Having obtained the recursion relations for the skew-orthonormal vectors, we introduce a convenient notation to 
express band-matrices P^^'' and R'^P'>: 



n-\-d 



n-\-d 



ra—n — d 



k=-d 



k=-d 



(2.12) 



and 



n-\~d 



n-\-d 



k=-d 



(2.13) 



where 

diag(^f (0), . . .,vf\n), . . .), := diag(Cr(0), . . . ,Cf ^(n), . . .), ^ diag(af (0), . . .,af{n), . . .). 

(2.14) 

Here, we have suppressed the x-dependence of a''^\n). This notation will be useful in deriving the folding function. 



More explicitly, with this notation, we can write 



and 



R 



Vn, 



Cf>(2) 


d^(0) 

cr(i) 

Cr(2) 






C!f\d+l) 








r^^(o) 


^i'^r^(O) 
^^(1) 




vT\id) 













x), the quaternions 


Tdin) . 






Cf (d + 2) 



d'',V2(2) rAi(2) Ci"j(2) 



>(/3) 








^r^d) 
^^(2) 



d^'jW ... 

v'-'kD ^id) ... 

d^r(2) d''j(2) d«(2) 



r?f (d + 2) 

We note that for even V{x), the quaternions of the outermost band are such that 
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l±d 



From (|2.6p . we get the following quadratic relation between the coefficients (Q , ?7j ) 

d~i 



d-l d~l 
■J=l 3=1 



^t'MK'i^ -I+J) -l^Cr,{n)vf\n -l+j) = C^W- 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



This is the compatibility relation for the difference-differential system and will be used in Section IX. 

Remark. Here, one can notice the basic difference between the properties of bi-orthogonal polynomials and their 
skew-orthogonal counterparts. In the skew-orthogonal vectors, the semi-infinite matrices are symmetric around the 
principal diagonal. This is not true for the bi-orthogonal polynomials where the matrices have an asymmetry. 
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3. THE DEFORMATION MATRIX 



In this section, we will consider infinitesimal deformation corresponding to changes in uk, the coefficient of the 
potential V{x). Using the definition of the semi-infinite matrices <i>^'^-'(a;) and $'^'^'(2;) (and ^''^'^^(a;), ^'^'^^(x)), the 
deformation matrix can be defined as: 



d 

duK 

d 
du 



K 



duK 
d 



du 



K 



K J ' 
D 

K J ' 



if = l,...2d, 
K =l,...2d. 



The matrix J/]^-* is anti-self dual, i.e. 



(/3) 
K 



(/5)l 



(3.1) 
(3.2) 

(3.3) 



Moreover, they satisfy the following relations with P and R (we drop the superscript /3 from the matrix Q^'^\ P*-'^' 
and R^^^ for simplicity): 



OukP = [Uf , P] , OukR = [Uf , R] . 



Explicitly, u'^'^ can be written as 



= ((Q")+ - (Q")") - ^((Q")o - (Q")"o), 



(3.4) 



(3.5) 



where Q^, denote the upper and lower triangular quaternion matrix while Qo is the diagonal quaternion. 
Proof: 

Differentiating Ea. p.l2p w.r.t. uk and using Eqs. (|3.ip and p.2p . we get Ea. p.3p . Eq. (|3.4p follows by interchanging 
the operators d/dx and xd/dx with d/duK- Finally to prove Eq. (l3.5p . we start with the (2 x 1) normalized quasi- 
polynomial 



$(f)(. 



1 



-^(-)n„(x), n e N, 



(3.6) 



where V{x) is defined in Ea. (|1.2p and 



n„(x) 



T^2n{x) 

7r2n+i(a;) 



(3.7) 



is the monic skew-orthogonal polynomial, defined in Eq. (|1.9p . To save cluttering, we have suppressed the (3 dependence 
of 7r„(a;). We have used = 52n+i- Differentiating with respect to the deformation parameter, we get 



duK'^i^\x) = L^aUK(5i^))$i«(x) + ^^^duKllnix) - ^<Pi^\x) 



Also differentiating the skew-normalization condition, we get 

d 



K 



duK 



duK{^i{'\x))'fi^\x) + ^i{'\x)duK'fi^\x) 



dx. 



1 



0^ + (Q^) 



D 



The diagonal elements have the form: 



2K 



n.n 



(3.8) 



(3.9) 



(3.10) 
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Thus, we get 



(3.11) 



Remark. Here the operation (Q^)^ means that the dual is taken first and then the lower triangular part collected. 
Also, due to the arbitrariness in the definition of an anti-self dual matrix, we may choose the lower off-diagonal 
element of the diagonal quaternions in C/]^'' zero. With this choice, we can remove the arbitrariness in the definition 
of the odd skew-orthonormal polynomials. 

4. THE DIFFERENCE EQUATION 

In this section, we introduce a sequence of companion- like matrices A'"j^\x) and B^^\x) of sizes {2d) x {2d). Using 
the relations in Eqs. (|2.4p and (|2.5p . we can write: 

d d 



k=-d 



and 



k=-d 

we get for N > d: 



a a 

E {ai!\n) - 6o,k) ^i%{x) = E {"k\n) - 6o,k) ^0 Vfc G N, 



and 



V '^N\d{^) j 



i ^Tdi-) \ 

V ^N+d-li^) J 

V 'i'N+d-li^) J 



where (suppressing the /3 dependence of aj{n)) 
( 1 











\-{a.i{N))-^Oii{N) ... -(a_d(A^))"'(«o(A^)) ... -(a-d(A^))"'«_d+i(A^) / 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



Remark. For potential V{x) with even degree and Vn, a±d{n) — R^J^\^^ — xP^^^^^ and its inverse has the quaternion 
structure: 



a±din) 



a 
b c 



Q:±d{n) ^ 



i 

a 

h_ 1 

ac c 



a, c 7^ 0. 



(4.6) 



Hence they are invertible. This is the criteria for the existence of the skew-orthogonal polynomials corresponding to 
even potential. We leave it to the reader to verify that a±d{n) is non- invertible for potentials with odd degree. 
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Using this, one can easily see that 



and 



(4.7) 



(AW,(a;))-i*W(:,), 



(4.8) 



where 



V 



Similarly, for ^'^(a;) and $(i)(x), A'^^\x) is r eplaced by i3('^'(x), where 
/ 1 

<(x) = 



V -(a_d(iV))-iarf(iV) ... -{a^d{N))-\aQ{N) ~ I) ... -{a-d{N))-^a^d+i{N) J 
More generally, this ladder operator can be used successively to obtain: 



(4 



3) 



(4.10) 



and 



^^^lAx) = (<l,(x))-i . . . (<l,(x))-i$(^)(a;), Vj e N, 



This is the underlying idea behind folding which will be discussed in details in the next section. 



(4.11) 



(4.12) 



FOLDING 



Remark. From this point on we focus on the skew-orthogonal vectors '^''n\x) (and its dual \E'n ^(x)), but everything 
being said can be immediately extended to the vectors (^n\x), (and hence ^'"n\x)) by interchanging the roles of 
matrices by R'^^'^ - 1 . 



The notion of "folding" is the following: we express any quasipolynomial <^\i^ [x) (or '^'^n' {x)) as a linear combination 
of 2d fixed consecutive vectors ^'^^ {x) = $^l_j(a^) , j = d, • ■ ■ — d + 1 with polynomial coefficients. We now provide a 
way of computing the folding function on the same line as done in [20l | for bi-orthogonal polynomials. 

For a fixed subsequence (or window) W , such that W ^ {N — d, . . . , N + d — 1},N > d, we seek to describe the 
folding of the infinite wave- vector ^^^\x) onto the window W by means of a single quaternion matrix F(x) of size 
oo X {2d) with polynomial entries such that: 



N+d~l 
k=N-d 



(5.1) 
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(Similarly for ^i^\x)). In fact it is more convenient to think of F(x) as a oo x oo matrix with only a vertical band 
of width 2d of nonzero entries (with column index in the range N — d, . . . , N + d — 1). 

Now, we will find an explicit form of the "folding function" in terms of some upper and lower triangular matrices. 
We will derive the folding function for ^'-^'>{x) and ^'(^^(a;). For the folding function for ^'^^^(2:) and $'^^^(0;) we 
follow the same procedure and hence state the result without proof. Note that everything will be written in terms of 
quaternion matrices. 

We will start with a few definitions and identities. We will introduce the shift matrix A comprising of 2 x 2 blocks 
of unit matrix: 



/O 1 ...\ 
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1 
1 



(5.2) 



satisfying the relation 



(5.3) 



i.e., in this case, the transpose is the same as dual. Hence for convenience, for the rest of the paper, we will denote it 
as A*. 

Let us define the projection matrices: 



N 



N 



n. n-n n 



(5.4) 



N 



M N AI 



We will also use the following identities. 



N+d N 



N+d 



N+d 



N+d 



N 



N+d 



N 



(5.5) 



N 



N+d 



nn-n n-n. (a*)V^i n. ^m-^- 



N+d N+d N N 



(5.6) 



To find an explicit formula for the matrix F{x) we will use the diagonal band-matrices aj, j = —d, ... d to express 

i?(/3) _^p(/3). 



k=l 



k=0 



(5.7) 



where A is the shift matrix defined in Eq. (|5.2p . (For simplicity, we have again suppressed the /? dependence of aj.) 
Using the shift-matrix on Ea. (|5.7p . we get 



d~2 



k=l 



k=0 



Thus we can define a semi-infinite matrix G^^\ which is strictly lower triangular: 

G(« ■.= l-l[-{Ara-_l[R''\x)]. 



(5.8) 



(5.9) 



Similarly 



A''a-^[R^'^\x)] = l + A'^a-i 



(5.10) 
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allows us to define a strictly upper triangular matrix: 



Multiplying Eq.dS^]) by a_d(A)'^ on the left and (1 - G^'^)) ^{K*faZ]i ™ the right, we get 
while for the upper triangular matrix, we have the relation 



(5.11) 



(5.12) 



(5.13) 



Thus the matrix [R'^^'^ — xP'^f^'^] has a left and right inverse, but they are not the same. One is upper triangular 
while the other is lower triangular. However, they satisfy the following relation: 



(5.14) 



One observes similar features for the matrix {Q — x) in the context of the two-matrix model studied in [20| . 

Following the above procedure, we can also get the upper and lower triangular matrix, needed to fold ^'^'^\x) and 
$(i)(a;). They are 



:= 1 - n -{^'fo.-_\[R^'\x) - 1], C(« := 1 - A'^a, " 1], 



such that they satisfy 



1 = [\r''\x) - 1]J [(1 - &^^r{K^fa-_\\ , 1 = (1 - C^^r AV[R''\-) - 1] 

A. Upper folding 

From the definition _ a;P(4)]$(4) ^ q, and Eqs. (EH) and we get 



(5.15) 



(5.16) 



Similarly from the definition [R'^^'> ~ xP^^'^\^^^\x) = 0, and Eqs. (fOj) and ([5?TT|) . we get 



\ d+lj 



(5.17) 



^^^\x) = &^^^^^^\x), ^^^\x) = Ig^^^^ + W\¥^^\x). (5.18) 

V d+lj 

Given that the quaternion matrix G is strictly lower triangular, we can use Ea. (|5.18p in an iterative manner to get 

N N N N 

n <f <'^(a^) = (1 - G(4^)"' n ^'"^ n *^'^(^)' n *^'^(^) = - ^^'^)'' n ^^'^ n *^'^(^)- (^.w) 



N 



N 



We will drop the parameter (3 for simplification. Also since the matrix G has d bands below the principal diagonal, 
we have for the folding function: 



N 



N 



\{F{x)^{l~G)-'\{G\{, 



(5.20) 



N 



13 



where N = N + d — 1. Simplifying further, we replace G from Eq. (|5.9p to get 

Iv Iv 



n^(:^)=(i-G)-^n 



Then using the identities obtained before, one gets 



N 



d 



N 



n 



(5.21) 



-I N 



i[F{x) = -(i-G)-^n"=^(^')'[^^''^(^)]n 



N 



N 



N 



N 



Thus we have for the upper folding 



N 



N 



N 



l[F{x) = -l[-{l-G)-'aZl{A') 



N 



B. Lower Folding 

Now, we will study the lower folding function. We can write 

M 



AT 



(5.22) 



M 



Y[<f^''Hx) = {1 - C^'^y) n^^'^n*^'^(^)' Y[^^'Hx) = il-C^'^f Y[C^'^Y[¥'\x). (5.23) 



M 



M 



M 



M 



where M = N — d. This relation is valid since we are folding finite band matrices. Thus F[x) can be written as 



M 



Y[Fix)^{i-cr'Y[cY[. 



(5.24) 



M 



Replacing C from Eq. (|5.1ip and using the identities, we get 



l[F{x) 



M 



M 



-il-C)-'A'a-^Y[R''\x)Y[, 

N 

N 

-l[+(l-C)-'A'a-'[R^''\x),Y[], 

M N 
N 

-l[+il-C)-'A'^a^'[R^'\x),Y[]. 



(5.25) 



M 



N 



C. Folding Function 



The folding function is 



JV 



N 



F{x)=l[Fix)+l[Fix)+l[, 



(5.26) 



M 



M 



14 



such that we get 



Fix) ^ ((1 - Cy'A'^a-' - (1 - G)-'aZliA'f) [R^''\x),Un]- 



(5.27) 



One can easily repeat the above calculations to find the folding function for 'I'('')(x) and <i>(^^(a;). It gives: 



Fix) ^ ( (1 - C)- A%-i - (1 - G)- aZliA')') [i?^^^(x),n^]. 



(5.28) 



Remark. In the context of bi-orthogonal polynomials, Bergere and Eynard [371] have found a very elegant expression 
for the folding function. It would be interesting to see the existence of such compact expression in the context of 
skew-orthogonal vectors. 



6. FOLDED DEFORMATION MATRICES 



In this section, we will show that a finite sub-sequence or window of these skew-orthonormal vectors satisfy a system 
of PDE's under the infinitesimal change of the deformation parameter. We will prove for $('^^(a;) and \E'*^^)(a;). It can 
be trivially extended to \l/('*)(a;) and <I>*^^)(a;). We state the result. 

For convenience, in this section, we will drop the superscript /3 from Q^^"^ as it has no importance on the deformation 
system. We define the folded deformation matrix: 

Iv Iv 

Y[ul^''Fix) ■.= U^''\ l[ul,'^Fi^) ■.^Uf'\ K^l,...2d, (6.1) 

M M 
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where we project a 2d x oo quaternion matrix Hm onto a oo x 2d quaternion matrix F{x) and F{x). We have 



M 
N r 



n 

M 



[Q?-(g^)e] + ^[Qo^-(Qo^)^ 



l[F{x) + f[F{x) + l[ 



N 



M 



M 



= - n l^^K^] n + n n - ihq'')'' n ^(^)' 

M MM MM 



M 
N 



M M 



= -nKin-n^" 



M 



D 



M 

Iv 



M M 

N 



N 



r''\x),1[ 



N 



.J|+(l_C)-i(a:X 

M 

N 77 Iv Iv 



AT 



AT TV 



nKin n^^n+nwvn-n 



M 

Iv 

n 

M 



MM N M 



{Q^''){l-C)-\a-_\K)\p^^KX{] 



M M 



{Q''){l-G)-\K'a-_\f 



r''\x),]1 



N 



N 



M 



N 



M JV+1 
1 



JV+l 

Iv 



N N 



n-n 

M N 

l[{Q^-x^f(R''\x))] 



r''\x),1[ 



JV+l 




N . 


M 






g (Qo^ + (Qo^)'') 

M 2 


n n 

M N 


iV 
M 


N . 


-2x^ 


TV N 

n-n 

M N 



N 



where 



5W, 



Wk{x) = {Q'' -x''){R''"{x))^ +{Q'' -x'^r{R''"{x))^ 
Now, if we define in analogy to the PauH matrix as, 

N N 



:= 



n-n 

M N 



= diag(l^,^l_^), 

d d 



(6.2) 



(6.3) 



(6.4) 



then we may write 







11m 




^^''^(^),n;v^ 


- 2x^E3, 


For ^^'^\x) and <^^^\x), we have 






11m 


(Qf + (5)°)5,3_n-VV,(.) 




- 2a;^S3, 



(6.5) 



(6.6) 
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where 



(6.7) 



Thus, we see that these skew-orthogonal vectors satisfy a systeni of PDE's with respect to the deformation parameter 
and is given by 



d 



duK 
d 



(6.8) 
(6.9) 



where $|^'(x) and ^'[^^(a;)is defined in Eqs. p.26|1.27p . 



7. THE DIFFERENTIAL EQUATION 

Finally, using the results of the deformation equation, we will derive the ODEs satisfied by the finite subsequence 
of skew-orthogonal vectors. 

Remark. We will calculate explicitly the differential operator for ^'^^^{x) and ^'^'^^(2:). The corresponding dual 
vectors ^'(^^(x) and <i>(^)(a;) respectively can be obtained by a simple transfer of Q*^^^ i-^- Q'^^^ and hence P'^' 1— > P^'^'' 
and i?'^^-' ^ R^^\ Everything else remains the same. 

We recall that the differential operator on any of these vectors: 

^$(^)(x) = P(«$(^')(x) = -[iV\Q))+ {V\Q)t + lliV'm, - iriQ))^]]<i>^PHx) 



2d-l 



K=0 



and similarly for ^'(''^(a;). Since 

d'}^^ and d}J^^ is the folded version of p(/5) , we observe that the differential equation 



(7.1) 



satisfied by the finite subsequence of $(^^(2:), namely $[^'(x), is given by 



M 



such that with Q := Q^^\ P := P^") and R := P^^)^ have 

N N 



M 

2d-l 



N 



2d-l 



Fix), 



K=0 



N 



2V'ix)E3 + Y[wf[ix) 

M 



R{x),\{ 

N . 



In the last step, we have used the anti-self dual property of P. 
Thus we have 



= Y{liPF{x)<^^'Kx)^ 2V'{x)^, + Y{lwUx)[R{x),Y{^] 



dx 



(7.2) 



(7.3) 



(7.4) 
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where 



Wi^(^) = [V'{Q)-V'{x)]{xP-R)-' + [[V'{Q))-V'{xr{xP-R) 



D , 



(7.5) 



Replacing <^'^^\x) ^ ^'^^\x) and Q^^) ^ g(i)^ (^nd hence pf^) 
differential equation for the dual vector ^^^\x). 
Explicitly, we have 



P(i) and i?'-*) ^ P(i) ) we can obtain the 



D 



(4) 



= 2t/'(a;) 



/l 













\ 

























1 



















1 

























\o 
















/ 















W^{x) 











2dx2d 








'ai(iV + l) 
-ad{N) 



a-d{M) \ 

: '•• 

a-i{N) ... a^4N) 





For vl'l^'' (x) we have 



UZpf{x)^^^H^)- 2V'ix)i:, + uZWKi^)[Ri^),UN] "fwi^) 



where }V^{x) is given by 



W^ix) = [V'iQ)-V'ix)mx)-l)j,: + [{V'{Q))-V'{x)r{R{x)~l) 



(7.6) 



(7.7) 



(7.i 



For the dual vector <^\J{x), we follow the same procedure. Replacing '^^^\x) ^ ^^^\x) and Q^"*) i— > Q^^-*, (and hence 
p(4) ^ p(i) and P('') P(i) ) we can obtain the ODE's for the dual vector $(i)(x). 



8. FUNDAMENTAL SOLUTIONS 

We have seen that a finite sub-sequence of (2 x 1) skew-orthogonal vectors $j^l^.(a;) and ^'j^l^(a;), fc = — (i+ 1, . . . d, 
satisfy a system of differential-difference-deformation equation. To obtain the fundamental solutions, useful in the 
Riemann-Hilbert analysis, we need to look for 2d other solutions. In this section, we will show that the Cauchy-like 
transforms of these skew-orthogonal vectors indeed form such solutions. 

We will look for some integral representations of these skew-orthogonal vectors as the possible "other solutions" . 
From the previous experience in orthogonal [10] and bi-orthogonal polynomials [llj, we look for some Cauchy-like 
transforms ^^\x) and l>^f '(a;), and hope that at least for sufficiently large ri, these functions will form the simultane- 
ous solutions for the system of difference-differential-deformation equations. The remaining solutions can be obtained 
by taking the Cauchy-like transforms of certain moment functions fj{x) of order j, to be defined below. This trick 
had been used previously to define Riemann-Hilbert problems for bi-orthogonal [ssj and skew-orthogonal polynomials 
[sot respectively. 

Let us define 



/, (x) :=exp[F(x)] / e{x-y)y'eM-y{y)]dy. J = 0, . . . , 2d - 2. 
Jk. 



(8.1) 
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We will show that the functions: 



*0-,/3)(a;) := /■ ZlI^L_Mfl^(f)(2)exp[-y(z)]dz, j = 0, . . . , 2d - 2, 

JR X Z 



(8.2) 

_ . . (8.3) 

X — z 

and the Cauchy-like transform of the skew-orthogonal vectors: 

¥fHx) = exp[T/(x)] / ^'^'\fM-Viz)] ^^^ ^ ^^^j^^^^j r ^i^\z)eM~Vi^)] ,^^ ^ ^ 

are simultaneous solutions of the differential-difference-deformation equations. 

A. The Cauchy-like transform 

Componentwise, we get using Eqs. (|1.19p and (|1.20p : 

i'i'^ (^) = E Pn%4>l^^ (^), 4>i'^ (^) = E (^)- (8-5) 

We also get 

;^^(«(x) = exp[F(.)] / - ^'(^) exp[-V-(z)]dz + E 



dx ^ X ~ z 

2d-2 

eMV{x)] E C(fHx)Z„,,„ +E^il^L^n2^), (8.6) 



m=0 



where y'(a;) = £V{x) and t/'(z) = j;V{z). Here, we have used eyip[~V{z)][{V' (x) ~ V {z)}/{x ^ z)] 

m— 



E'to C^ri\x)4>\^\z)- Thus we have 



^^i«(x)=E^i«H^), n>2d. (8.7) 

Using Eqs. ((83)) and ((87)) . we get 

'^f^(^)=^i'^(^)=E^"!™'^™ (^)' V;f^(:^) = 20«(x), 0:«(a;)=E^i:i'^^H^): «>2d. (8.8) 

Thus for n > 2d — 1, our chosen functions satisfy the same system of differential equations as done by the skew- 
orthogonal vectors (j)]^\x) and ipi!^\x) (|2.ip . 

Having obtained the relations between the 0'^'''(x) and 'ij)'^^\x), we will now show that these functions also satisfy 
the same recursion relations. We start with 

^Ul^H^) - 4exp[^(.)] / ^^fWIzM,,, 
dx dx 7jj (x — z) 

= xV'{x)^'i^\x)-eMy{x)] I ^^5^Vi'')(z)exp[-F(z)]dz + exp[T/(x)] / ^^^^^ exp[-F(z)]dz. 



dz 



X ~ z 

Now, replacing x^x — z + zon the right hand side, we get 
d 7,,),^, PT./^M [ ^V'{x)~zV'{z) 



x^¥n^\x) = exp[y(a;)] / ^^iA^L^^^(« (z) exp[-^(z)]dz + exp[y(x)] / (^(«(z) exp[-V^(z)]) 
ax .ha X — z Jo az \ / 



+ exp[l/(a;)] / 

JR 



exphF(z)] 



dz 



dz. (8.9) 
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Thus we have 



, 2d-l 

a 



X-^¥n^\x) = exp[^(x)] ^ SiXx)Z.^,n + ^ ^i^i^i'^H^) ' = 1' 



da: 

m— 
2(i-l 



= exp[F(x)] ^ d(^)(a;)Z™,„ + - 5r,^mWn^\x), /3 - 4. (8.10) 

Here, we have used exp[— F(z)] [{xy(a;) — zV'{z)}/{x — z)] = Y^mZo dm\x)(j)m\z) . Thus for n > 2d, we can see 

that the first term drops out and hence ^p'if ^x) satisfy Eas. (|2.21 12. 3p respectively. 
Similarly, for (3 = 4, we get for n > 2d, 

ax 



= exp[y(x)] / I^^±il^W(^)exp[-y(z)]dz, 

JR X ~ Z 

where the first term drops off due to the condition n > 2d, d > 1. For = 1, we differentiate (IS.lOp to get 



x^4>i'Hx) ^ Y.(^nl-Sn,m)4>i'Hx), n > 2d. (8.12) 



dx 

This in turn proves that for n > 2d, the Cauchy-like transforms satisfy the same set of recursion relations (|4.ip and 
(|4.2p as done by the skew-orthogonal vectors. 

Thus we prove that our function is compatible under a system of difference-differential equation, i.e. the solution 
of the system of differential equation is also a simultaneous solution of the system of difference equation. We will now 
show that these solutions are also compatible with the system of deformation relations i.e. they are also solutions of 
the system of PDE's. 

We can see that 

d ~(^^ , , exp[V{x)] f x^ -z^ 



^(«(.) . '-^^^ / ^^^^^i^Hz)eM~V{z)]dz + J2iU^^)_^l^Hx) 

(JUk J\ Jb X — z 



k-l 



eMV{x)]J2a(i\x)Z^,^ + J2{U^^)^^J(£\x). (8.13) 



m— 



where exp[-F(z)] [{x^ - z^}/(x - z)] = Y.fn=l a'i\x)4fj (z). 

Thus for n > 2d, 'ipl'!^\x) satisfy Eq. (|3.2|) . Also, combining Eqs. (|8.8| I8.13p . we can see that 



4>i^Hx) = E(^^)«,™<^™ (^)' "^2d. (8.14) 
ons 

equation and hence are mutually compatible 



Thus for n > 2d, the new-functions 4>n^\x) and ^pil^^x) satisfy the system of differential-difference-deformation 



B. The other solutions 

In this subsection, we will show that '^n''^\x) and ^n'^\x), j = 0, . . . ,2d — 2, are also solutions to the system of 
difference-deformation-differential equations. For this, we will characterize the skew-orthogonal functions (l)^i}\x) and 

V'i^'' (x) (which are polynomials of degree k and k + 2d — 1 respectively) through a set of skew-orthogonality relations 
with respect to the functions fj{x). The expressions for these auxiliary solutions were suggested (for the case /3 = 1) 
by the common referee of this paper and |39| using expressions contained therein, and extended by the present author 
to the case /3 = 4. 

We will start with the identity for a polynomial 7rj(x) of order j as: 
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TT,{x) = exp[T/(a;)]^(a;^-2rf+i exp[-y(x)]) (8.15) 

There are 2fc + 1 orthogonality property for the skew-orthogonal polynomials of order 2k and 2/c + 1. Starting with 

Zj+2d-i^2k ^ J^yj!ff}TT,+2d-iexp[-V{x)]dx = 0, Vj = 0,...,2fc-2d+l, (8.16) 
and using Eq. (|8.15p . we get the 2fc — 2c? + 2 conditions for 

/3 = 1, Z2kj+2d-i = 2 / (bi'ki^W exp[-Vix)]dx = 0, Vj = 0, . . . , 2fc - 2d + 1, (8.17) 

and 

/3 = 4, Z2kj+2d-i= [ i''2k\^>'('M-yix)]dx = 0, Vj = 0,...,2fc-2d+l, (8.18) 
while for the remaining 2d — 1 conditions, one can have 

Z2k.j^ f f ^^^^{x)y^eix-y)eM~Viy)]dxdy ^ 0, j = 0, . . . , 2d - 2, (8.19) 

Zj,2k^ I I i^^^^{x)y^+^'-\{x-y)eM-V{y)]dxdy = 0, j = 0, . . . , 2d - 2. (8.20) 

Thus we define the function fj{x): 

f,{x) := eMV{x)] f e{x - y)y' eM-V{v)]dy, j = 0, . . . , 2d - 2. (8.21) 

Remark: Here, we note that for /? = 4, without loss of generality, fj{x) can also be defined as: 

f,{x) := eMVix)] J f eM-V{y)]dy + c, j = 0, . . . , 2d - 2, (8.22) 

where c is a constant that can be fixed by skew-orthonormalization condition. 
Having defined fj{x), we will show that the functions: 

/ l2M^IM^i0)^,)eM-V{z)]dz; <I>(f''')(x) := / ZlM^^M $(«(^) exp[-y(z)]dz, (8.23) 



satisfy the same system of difference-differential-deformation equations. 
Componentwise, using Eqs. p.lOp and (|1.20p . we get 

We also get 
dx 



rn'-'^Kx) = f,{x) / ^Vr(^)exp[-^(z)]dz + 2 / ^W(^)exp[-F(z)]dz 



X — z 



ili^LJM^f (^) exp[-F(2)]dz 

X — z 

2d-2 j-1 



7n=0 m—0 
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Thus we have 



^i^i^'^Hx) =Y.Pi%i^l^^\x), n > 2d. (8.26) 



dx 

Using Eqs.dHU) and (E^ll), we get 

q^:['^'Hx) = ^(^■■^)(x) =Y.Pi%^ii^'Hx), ^i'^'\x) = 2</)(f^i)(x), 0;(^-i)(x) = E^n:i</'^''^(^), ^ > ^27) 

Thus for n > 2d, our chosen functions satisfy the same system of differential equations as done by the skew-orthogonal 
vectors (t>i^\x) and tpi^\x) ((TT|) . 

Having obtained the relations between the 0^^''^^ (x) and ijj'^^'^'> (x), we will now show that these functions also satisfy 
the same recursion relations. We start with 

x^^i=''\x) = f,ix) f ^""^'^f " ^"^V (^)(z) eM-~V{z)]dz + 2 / ^"'^ ' V i^^ (z) exp[-y(z)]d. 



(^(«(z)(/,(x) -/,(z))exphy(z)]) dz+ / zi^f\z y'^''>^ _ ^J^'> exp[-V-(z)]d#.28) 



which for /3 = 1 gives 



, 2d-l j 

^^^n-'H^) - /,(x) ^ dW(a;)Z™,„ + 2 5]a«(x)Z„,„+^i?lM^"-i)(x), j = 0, . . . , 2d - 2, (8.29) 

m— m— 

and for P — 4 gives 

, 2c;- 1 i 

a^^V'i-'''^ (X) = /,(^) E (^)^-." + 2 E (^)^™." + J2(^nl - 6n,n.)^^r^-'^^ (x) , J = 0, . . . ,2d ^8.30) 

rn—0 m—0 

Thus for n > 2d, we can see that the first two terms drop out and hence ipn'^\x), like tpr!^\x) satisfy Eqs. ()2.2p and 
(12. 3p respectively. 

Similarly, for /3 = 4, we get for n > 2d, 

= E^"1<^™'^W' (8-31) 
where the first term drops off due to the condition n > 2d, d > 1. For P — 1, we differentiate (|8.10p to get 

X^<t>li''\x) ^ Y.(^nl-Sn.nMl-'\x), ^ > 2d. (8.32) 

This in turn proves that for n > 2d, these functions satisfy the same set of recursion relations (|4.1I) and (|4.2p as done 
by the skew-orthogonal vectors. 

Thus we prove that our function is compatible under a system of difference-differential equation, i.e. the solution 
of the system of differential equation is also a simultaneous solution of the system of difference equation. We will now 
show that these solutions are also compatible with the system of deformation relations i.e. they are also solutions of 
the system of PDE's. 

We can see that 

^V'P^(-) - f ^-—^^i'H^^)eM-V{z)]dz + Y^iU^^l^i,ii^PHx), 



m—0 

Thus for n > 2d, tpi'''^\x) satisfy Eq. (|^ . Also, combining Eqs. (|8.27l8.33p . we can see that 

'l^'ii-^Hx) = E(f^^U™'^™''^(^)' "^2d. (8.34) 
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Remark: Here we encounter terms of the form fj+K{x), j = 0, . . . ,2d — 2 and K = I, . . . , 2d. For (j + K) > 2d — 2, 
this term does not contribute, by definition. Also if we consider only K = 2d, i.e. the leading term in V{x), no such 
problem arise. 

Thus for n > 2d, the new-functions 
equation and hence are mutually compatible 

Explicitly, we have 



n'^^ (x) and ipn'^'^ (x) satisfy the system of differential-difference-deforniation 



/ 



[W] 



*iv'L(^) *r.(^) 



>(/3) 



(0,(3)1 



(2d-2,/3) 
N-d 



V "^'nU-A^) 



(0,/3) 



(2d-2,/3) 



(/3) 



(/3) 



(x) <P^^U^) 



(0,/3)/ 



(2rf-2,/3) 
N-d 



(^) y 



(/3) 



,(0,/3) 



.35) 



Thus we have obtained the fundamental solution to the system of differential-difference-deformation equation, 
satisfied by the skew-orthogonal vectors. From here, it should be possible to define the RHP for these skew-orthogonal 
vectors. It has been obtained in Ref.[3§| for /3 = 1 and can be extended to /3 = 4. 



9. COMPATIBILITY OF THE FINITE DIFFERENCE-DIFFERENTIAL-DEFORMATION SYSTEMS 

In the previous section, we have obtained the Fundamental system of solutions for the differential-difference- 
deformation equation. 

In this section, we will provide an alternate proof to show that the recursion relations, the linear differential equations 
and the deformation equations are compatible in the sense that they admit a basis of simultaneous solutions, provided 
the vectors form a certain algebra defined by the commutation relations (|2.6p and (|3.4p . In other words, here, we will 
not use the explicit form of the fundamental solution to prove the compatibility. 

This is the same as saying that the shifts W W + 1 in Eqs. (|4.4[ ), (|4.8p respectively, are compatible as vector 
differential-difference systems. This means that there exists a sequence of fundamental matrix solutions and 

,(/3) 



'^[w](^) "^hich are simultaneous solutions of 



*[w+i]W = (9.1) 

^*[^](x) = <^^(x)*[^](x), (9.2) 

and 

*[w+i](^) = B^^\x)^^^^{x), (9.3) 

^*[^](x) = D^^\x)^[^^{x), (9.4) 

respectively. The same result holds ior P — 1. This means that there exists a sequence of fundamental matrix solutions 
satisfying the recursion and differential equations 

*[w+i](^) = a!^^(^)*[^](^), (9.5) 

^*[^](x) = i^«(x)*[^,(x), (9.6) 

and 

*[w+i](^) = B'-^\x)^[^^ix), (9.7) 

^*[^](^) = D}^\x)^[^^ix), (9.8) 

respectively. 

One must note that the differential equations are also compatible with the shift W ^ W — 1. The proof is exactly 
similar to the above case. 

Similar procedure will be repeated to show the compatibility between the difference-deformation equation. 
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A. Proof: Difference-differential equation. 



Now we will prove the compatibility of the differential-difference equation for ^^^^ (x). The others are exactly similar 
and hence not repeated. The proof follows on similar line to the one outlined for bi-orthogonal polynomials [l^ . 
Let us define a finite subsequence (or window) containing functions: 



which is a solution to the differential equation 



N> d. 



(9.9) 



(9.10) 



7ic(4) 



A glance through the window will confirm that the shift W i — > Vl^ + 1 will introduce one new function ^j^_^_^{x) 
inside the shifted window. Componentwise, let us suppose that they satisfy the recursion relation 



K4) 



[q^dim) ~ xC-dim)]<S?'Z+dix) ^ X 0(™)*m-/(a;) - m > d. (9.11) 

i=-d+l l=-d+l 



It is easy to see that the differential equation componentwise reads: 

dx 



^<!>i^\x)^ Y G(n)<'!,(x), n^N~d,...N + d-l, 



(9.12) 



where the (x) outside the window can be expressed in terms of that within the window by using Eq. (|9.11|) 
recursively. To show compatibility between the shift and differential operator, we need to show that the newly defined 
function, which is assumed to satisfy 



:5c(4) 



(9.13) 



:f:(4) 



will also satisfy the same differential equation i.e. we must show that the newly defined function $jy+(j(^) "^il^ ^-l^o 
be a solution to the differential equation: 



l+d{ 



(9.14) 



Remark. One must be careful about the position of the inverse of the quaternion, as in general, they do not commute. 

(4) 

Once this is proved, we can argue by induction that ^N+d+ji^) satisfy the same sort of differential equation for 
any j > 1. To do this, we compute 

'^-'^^^^(^)^^'+'^('^^ ^ ^ ry,(7V)¥iv^(:r)') 



l — — dj — — d 
r(4) 



i=-d 

d d 

-EE m{N)UN~m'^ii^,{x 

l— — d j—^d 
d 

+V-d{N) E Ci{N + dW^^li+di^). 

l = -d 

From Ea. (|9.15p . we can see that in order to prove Eq. (|9.14p . we will have to show that 



(9.15) 



E ci^nU^) + E E o(A^)^.(A^ - mTi^^i^) -Y. J2 viiNXA^ - inTi-A^) - o. (9.16) 



^:f:(4) 



:f:(4) 
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But this is nothing but a direct consequence of Eq. (|2.18p . in terms of its components. We can repeat the same 

(4') 

procedure for the shift W i — > W 1. FinaUy, one can extend this argument by induction to <I>jy^^(x), thereby 
completing the proof. 

We can repeat the same procedure to the sequences 



■tt7(1)' 



Tt7(1)' 



*W-d(2;),--->*Ar+d-l(a;) 



N>d, 



(9.17) 



and with some minor modifications to 



7jr(l)/ ^ r7r(l)* / x 



(1)* 

,$Ar+d-l(a;) 



N > d. 



(9.18) 



and 



N>d, 



(9.19) 



to prove compatibihty conditions between the difference-differential equations satisfied by the skew-orthonormal vec- 
tors. 



B. Proof: Difference-deformation equation. 



We show that these functions are also compatible with the deformation equations. We will briefly outline the proof. 
We start with the same finite subsequence (or window) containing functions: 



$l^(a;):= <^f^_^{x),...,<S>N+d-iix) 



N> d. 



which componentwise, will satisfy the PDE: 

K 



duK 



n^N~d,...,N + d-l. 



Then under a shift W i — > + 1, we get the new function ^N+di^) defined by 



which should satisfy 



d -(4) 
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l=-d+l 



duK 



^nU^) = E + d)^V+d-,{x)- 



To do this, we compute 
/ d 
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duK 
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J2 Ci{N)<p'^Ux) ~ E m{N)^NUx) 
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d K 

l = -d l=-dj=0 

E ^iiiN)^N\{x)-il Y.m{N)U^{N~m''^_^_,{x) 
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3=0 



(9.20) 
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(9.24) 
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where ' denotes . Rearranging the coefRcients in Eq (|9.14p , we can see that basically we will have to show that 

d d K 

i=-d i=-dj=a 

d d K 

- ^[{N)^%Ux) - Y.^muf{N ~ (9.25) 

j=-d l=-dj=0 

But this is nothing but a direct consequence of the string Eq. ()3.4|) . in terms of its components. Hence we get Eq. ()9.23|l . 
We can repeat the same procedure for the shift W i — > W — 1. Finally, one can extend this argument by induction 

(4) 

to thereby completing the proof. 

We can repeat the same procedure to the sequences ^fj^^x) and with some minor modifications to (x) 

(4) 

and {x) to prove compatibility conditions between the difference-deformation equations satisfied by the skew- 
orthonormal vectors. 

10. CONCLUSION 

In conclusion, we have obtained a system of differential-difference-deformation equation for a finite subsequence 
of skew-orthonormal vectors. We see that similar to the orthogonal and bi-orthogonal polynomials in ordinary 
space, the Cauchy-like transforms of these skew-orthogonal vectors, of sufficiently high order (to be precise, N > 
2d) also satisfy the same differential-difference-deformation equations in the quaternion space. We also derive an 
integral representation of these skew-orthogonal vectors in order to obtain the fundamental system of solutions for 
the overdetermined system of ODEs, difference and deformation equations. 

On the other hand, from our already existing knowledge of skew-orthogonal polynomials, we know that the final 
results of the RH-analysis will have some distinct differences from the orthogonal polynomials. For example, the zeros 
of orthogonal polynomials are real, which is in contrast to that of the skew-orthogonal polynomials To understand 
these properties, one needs to extend the already existing theory for orthogonal polynomials and matrix-RHP [IJJ to 
skew-orthogonal polynomials and q-matrix RHP. We wish to come back to this in a later publication. 

Finally, we would like to mention a result of Dyson and Mehta. They showed (49| that for the circular ensembles, 
"the probability distribution of a set of N alternate eigenvalues of a matrix in the Orthogonal ensemble of order 2N is 
identical with the probability distribution of the set of all eigenvalues of a matrix in Symplectic ensemble of order A^" . 
It would be interesting to know if something similar exist for ensembles with polynomial potential and the role (if any) 
played by the duality relations between different skew-orthogonal vectors of Orthogonal and Symplectic ensembles. 
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